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We describe calculations of bb production to next-to-next-to-leading order (NNLO) and 
next-to-next-to-leading logarithm (NNLL) near threshold in pp interactions. Our calcu- 
lations are in good agreement with the bb total cross section measured by HERA-B. 
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Oh: 

Factorization assumes it is possible to separate QCD cross sections into universal, non- 
perturbative parton densities and a perturbatively calculable hard scattering function, the 
partonic cross section. However, some remnants of long-distance dynamics in the hard 
scattering function can dominate corrections at higher orders near production threshold. 
These Sudakov corrections have the form of distributions singular at partonic thresh- 
old. Threshold resummation techniques organize the singular distributions to all orders, 
extending the reach of QCD into near threshold production. The singular functions or- 
ganized by resummmation are the plus distributions, [In x/x] + , where x denotes the 
'distance' from partonic threshold. 

The first attempts to resum the heavy quark cross section were at leading logarithm 
(LL) and exploited the fact that to LL, the Sudakov corrections to the heavy quark 
cross section were identical to those obtained for Drell-Yan production This early 

resummation calculation, like some of the later results that followed, used an empirical 
cutoff to keep the strong coupling constant from blowing up. Resummation beyond LL 
cannot make use of this universality because the color structure of each partonic process 
must be treated separately ||. The NLL QQ terms were first resummed in Ref. [|] for 
a simplified case and later fully solved for the qq channel 0. The exponents in the 
gg channel were calculationally unwieldly, requiring large cutoffs even for it production 
where the resummation should work best. However, one advantage of the resummed cross 
section is that when it is expanded in powers of a s , it provides estimates of unknown finite- 
order corrections without resorting to a cutoff or other prescriptions. We have calculated 
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the double-differential heavy quark hadroproduction cross sections up to next-to-next-to- 
leading order (NNLO), O(af), and next-to- next-to-leading logarithm (NNLL), i.e. keeping 
powers of the singular functions as low as I = 2i — 1 at order 0(al +3 ) where i = 0,1, . . . 
[0,|J. Since resummation is based on expansion of the LO cross section, we only discuss 
QQ production in the ij = qq and gg channels since qg scattering first appears at NLO. 

In our calculations, the distance from partonic threshold in the singular functions de- 
pends on how the cross section is calculated, either by integrating away the momentum 
of the unobserved heavy quark or antiquark and determining the one-particle inclusive 
(1PI) cross section for the detected quark, or by treating the Q and Q as a pair in the 
integration, pair invariant mass (PIM) kinematics. In 1PI kinematics, 

p(P 1 ) + p(P 2 ) — Q( Pl ) + X(p x ) , (1) 

where X denotes any hadronic final state containing the heavy antiquark and Q{pi) is 
the identified heavy quark. The reaction in Eq. (H) is dominated by the partonic reaction 

i(h)+j(k 2 ) — Q( Pl ) + xm(p> 2 ) . (2) 

At LO or if X[Q](p' 2 ) = Q(p 2 ), the reaction is at partonic threshold with Q momentum 
p 2 . At threshold the heavy quarks are not necessarily produced at rest but with equal 
and opposite momentum. The partonic Mandelstam invariants are 

s=(k 1 + k 2 ) 2 , U = (k 2 - pi) 2 - m 2 , u = (k\ - pi) 2 - m 2 , s 4 = s + U + u\ (3) 

where the last, = (p 2 ) 2 — m 2 , is the inelasticity of the partonic reaction. At threshold, 
S4 = 0. Thus the distance from threshold in 1PI kinematics is x = s^/m 2 . In 1PI 
kinematics, the cross sections are functions of t\ and u\. In PIM kinematics the pair is 
treated as a unit so that, on the partonic level, we have 

i(h)+j(k 2 ) — QQ(p')+X(k'). (4) 

The square of the heavy quark pair mass is p' 2 = M 2 . At partonic threshold, X(k') = 0, 
the three Mandelstam invariants are 

M 2 M 2 
s = M 2 , tx = — — cos0), ui = — — (l + £ M cos0) (5) 

where Pm = yl — Am 2 /M 2 and 9 is the scattering angle in the parton center of mass 
frame. Now the distance from threshold is x = 1 — M 2 /s = 1 — z where z = 1 at 
threshold. In PIM kinematics the cross sections are functions of M 2 and cos 9. 

The resummation is done in moment space by making a Laplace transformation with 
respect to x, the distance from threshold. Then the singular functions become linear 
combinations of ln fc N with k < I + 1 and iV = Ne lE where 7^ is the Euler constant. The 
1PI resummed double differential partonic cross section in moment space is 
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To find the PIM result, transform t% and U\ to M 2 and cos# using Eq. (|5|). (P) P refer 
to (anti-)path ordering in scale /i'. The cross section depends on the 'hard', Hy, and 
'soft', Sij, hermitian matrices. The 'hard' part contains no singular functions. The 'soft' 
component contains the singular functions and, from its renormalization group equation, 
the soft anomalous dimension matrix T l g, 2 dimensional for qq and 3 for gg, can be 
derived. The universal Sudakov factors, the same in 1PI and PIM, are in the exponents 
Ei, expanded as 



exp(^(iV u , / i,m))^l + ^f5:^ 1) ln fc (iV M )j + ^j [£ C£ (2) \n k (N u ) j + . . . . (7) 

The coefficients C% , as well as the detailed derivation of the resummed and finite- 
order cross sections, can be found in Ref. @. The momentum space cross sections to 
NNLO-NNLL are obtained by gathering terms at O{o i s ) and 0(af), inverting the Laplace 
transformation and matching the iV-independent terms in Hij and SV,- to exact 0(af) 
results. 

We have studied the partonic and hadronic total cross sections of tt and bb production. 
Any difference in the integrated cross sections due to kinematics choice arises from the 
ambiguity of the estimates. At leading order, no additional soft partons are produced 
and the threshold condition is exact. Therefore, there is no difference between the total 
cross sections in the two kinematic schemes. However, beyond LO and threshold there is 
a difference. To simplify the argument, the total partonic cross section may be expressed 
in terms of dimensionless scaling functions that depend only on rj = s/Am 2 - 1 ||, 
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a ij[S ,m-,^) = (W/,)) fe £ fS'Hri) In' f ^ ) . (8) 

" b k=o i=o x" 1 



We have constructed LL, NLL, and NNLL approximations to / i7 - ' in the qq and gg 



channels for k < 2, I < k. Exact results are known for k = 1 and can be derived using 
renormalization group methods for k = 2, I = 1, 2. Thus the best NNLO estimate of the 
cross section includes the exact scaling functions and the NNLL estimate of fjf'°^. On 



the left-hand side of Fig. [I], we compare fL' and fig'°\ the only approximate scaling 
function, in 1PI and PIM. The results are quite similar at small t] but begin to differ for 
i] > 0.1, especially in the gg channel. If the parton flux is maximized for r\ < 1, as for 
the HERA-B energy, y/S = 41.6 GeV, the reaction is close enough to threshold for the 
results to be reliable. Unfortunately at RHIC, the flux peaks at r) ~ 1, making predictions 
at RHIC from the threshold approximation alone unreliable. The gg channel dominates 
bb production in pp interactions. An inspection of the scaling functions shows that the 
results could differ substantially between the two kinematics. 

The total hadronic cross section is obtained by convoluting the total partonic cross 
section with the parton densities ff evaluated at momentum fraction x and scale /i, 

a pp (S,m 2 )= f dr [ —ff(x 1 ,iJ, 2 )ff(—,fxAa ij (TS,m 2 ,iJ 2 ). (9) 

■ • - J4m 2 /S Jt X\ \X\ J 



If the peak of the convolution of the parton densities is at rj < 1, the approximation 
should hold. On the right-hand side of Fig. |l] we compare the 1PI and PIM results with 
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Figure 1. The left-hand side shows the ^-dependence of the NNLL scaling functions. 
We show ffffi°\r)) in 1PI (solid) and PIM (dashed) kinematics and ffjg ,0 \r)) in 1PI (dot- 
dashed) and PIM (dotted) kinematics. The right-hand side compares the total bb cross 
sections at fixed target energies calculated with CTEQ5M and fi = m = 4.75 GeV. The 
exact NLO result is shown in the solid curve while the NNLO-NNLL results for 1PI and 
PIM kinematics are given by the dashed and dot-dashed curves respectively. 



the exact NLO results, all calculated with the CTEQ5M parton densities || and \i — m. 
The NNLO-NNLL corrections are substantial. The average of the NNLO-NNLL 1PI and 
PIM cross sections, 

(7^(41.6 GeV) = 30±8± 10 nb , (10) 

is in good agreement with the bb total cross section measured by HERA-B |10|| . The first 
uncertainty is due to the kinematics choice, the second to the scale dependence. The 
uncertainties in scale and kinematics choice are essentially equivalent. 
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